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We investigate the potential energy surface of a (j> 4 model with infinite range interactions. All 
stationary points can be uniquely characterized by three real numbers «+, qq, Q— with q+ + qo + 
a_ = 1, provided that the interaction strength u is smaller than a critical value. The saddle index 
n s is equal to ao and its distribution function has a maximum at n™ ax = 1/3. The density p(e) of 
stationary points with energy per particle e, as well as the Euler characteristic x( e )i are singular at 
a critical energy e c (/x), if the external field H is zero. However, e c (pi) ^ v c {(j), where v c (fj.) is the 
mean potential energy per particle at the thermodynamic phase transition point T c . This proves 
that previous claims that the topological and thermodynamic transition points coincide is not valid, 
in general. Both types of singularities disappear for H ^ 0. The average saddle index n s as function 
of e decreases monotonically with e and vanishes at the ground state energy, only. In contrast, the 
saddle index n s as function of the average energy e(n s ) is given by n s (e) = 1 + 4e (for H — 0) that 
vanishes at e = —1/4 > Vq, the ground state energy. 

PACS numbers: 05.70.Fh, 61.20.Gy, 64.70.Pf 
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I. INTRODUCTION 

Topological features play an important role in sev- 
eral branches of physics. Examples in condensed mat- 
ter physics are discussed in Ref. 0. Those examples 
do not include thermodynamics and phase transitions. 
That topological concepts might be relevant for equilib- 
rium phase transitions has already been emphasized long 
time ago @ and that they can be very useful in condensed 
matter physics has been demonstrated recently 0, 0] . 

Usually, equilibrium phase transitions are indicated by 
a singularity at the transition temperature T c of thermo- 
dynamic quantities, like free energy, specific heat, etc. 
One may ask whether other indications for such phase 
transitions really exist. This question has been studied 
by several groups in recent years. Geometrical entities 
like the Ricci curvature, and dynamical ones like Lya- 
punov exponents were used for a classical planar Heisen- 
berg model with nearest neighbor interactions and di- 
mension d = 2,3 a nearest-neighbor 4 model for 
d = 3 with 0(n)-symmetry (n = 1, 2, 4) 6] and in d = 1,2 
with 0(l)-symmetry 0- It has been found that both 
quantities exhibit a singularity at a critical energy per 
particle, e c , for those dimensions for which an equilib- 
rium phase transition occurs at T c > 0. Furthermore, 
v c , the internal energy per particle at T c , equals e c , i.e., 
the geometrical and thermodynamic singularityoccur at 
the same energy. It has also been speculated 0, Q that 
these singularities are related to qualitative changes in 
the topology of the potential energy surface (PES) of 
those models. That this is true indeed has been proven 
first for a mean-field XY-model @- 

One of the most interesting topological quantities is 



the Euler characteristic X S which is a topologi- 
cal invariant. For the two-dimensional nearest neighbor 
</> 4 model [l(J, the mean-field XY and mean-field k- 
trigonometric model [l2l[T^ | it was proven that x(e) also 
becomes singular at e c = v c . In addition, it was shown 
[T2I IT^ | that the type of singularity depends on the order 
of the phase transition. 

X(e) is directly related to M(e, N s ), the number of sta- 
tionary points of the PES of a A^-particle system with 
energy < e and saddle index n s — N s /N^, 



N 



X (e)= £(-l)".M( e ,JV a ). 



(1) 



N„=0 



N s is the Morse index, i.e., the number of negative eigen- 
values of the corresponding Hessian matrix of a station- 
ary point. M{e,N s ) is an exponentially large number in 
JV» 1: 



M(e,N s ) ~exp[Ns(e,n s )}, 



(2) 



where s(e, n s ) is the configurational entropy per particle. 
The corresponding density of states p(e,n s ) is given by 



P(e,n s ) 



d_ 

de 



M{e,N s ) - exp[Ns(e,n s )}. 



(3) 



The relationship between Eqs. and leads to the 
assumption that the singularity of x( e ) may originate 
from a specific behavior of p(e, n s ) for N — ► 00 or that of 
s(e,n s ). 

The density of states p(e,n s ) plays an important role 
in the investigation of the PES. For instance, one can 



2 



define the average saddle index for a fixed energy 

n s (e) = / dn s n s p(e,n s ) (4) 
Jo 

and the average energy for a fixed value of n s 

/oo 
de e p(e, n s ). (5) 

In the limit N — ► oo the averages n s (e) and e(n s ) are sim- 
ply the values that maximize the configurational entropy, 
i.e., they are solutions of the equations 



ds . 

-^—{e,n s 
un s 



= 0, 



|(e,n s ) = 0. 



(6) 



The saddle index properties of a PES have also played 
an important role in another respect. Studyin g gl assy 
dynamics and ideal dynamical glass transition f 1 - ll ] for 
liquids it has been found numerically [TH Il6j| that the 
temperature-dependent average saddle index practically 
vanishes at a temperature T* , which is close to the mode- 
coupling glass transition temperature T c jljj]. However, 
this conclusion should be taken with some care. First 
of all most saddle poin ts were quasi saddles (see the 
discussion in Refs. |l7l IlsL IT^ I and second, plotting 
log n s (T) (or a related quantity) versus \JT does not 
exhibit a quasi sing ular behavior at T c [20, ij. That 
n s (T) vanishes at T = 0, only, has been proven for the 
/c-trigonometric model [13| . This indicates that n s (T) for 
systems with self-generated disorder behaves differently 
than for systems where the disorder is quenched. For 
the latter it has been proven that fi s (T) vanishes at the 
dynamical transition point, at least for mean field like 
models [22|. 

A numerical determination of the true saddles of a bi- 
nary liquid with particle number N < 13 gives evidence 
that fjft(e) vanishes at an energy e*, which still depends 
on N jlj]]- This evidence holds for both cases where the 
average n s is plotted versus e and n s is plotted versus e 
|17|. Whether this vanishing at e* is spurious or not is 
not known. It is obvious why vanishing of n s at T* or 
e* may be relevant. In that case, the system is mostly 
close to local minima (n s = 0) for T < T* or e < e* 
and the dynamics is dominated by activated processes, 
in contrast to T > T* or e > e*, where the particles dy- 
namics is more flow-like. Hence, vanishing of n s (e) may 
indicate a qualitative change in the dynamics fijj • 

There is another result presented in Ref. |l7j| which 
concerns the distribution function of the saddle index 
P(n s ): 



p(n s ) = / dep(e,n s ). 



(7) 



It is found that p{n s ) is a Gaussian with a maximum 
at n™ ax w 1/3. Although the physical relevance of this 



result is not clear, it seems to be an interesting property 
of the topology of the PES. 

We hope that the exposition above has made obvious 
the role of topological features for both the thermody- 
namic and dynamic behavior. It is the main purpose of 
our paper to analytically investigate for a mean-field 4> 4 
model the existence of a singularity in the topology of 
its PES and the relation to a thermodynamic singular- 
ity and to calculate the saddle index properties discussed 
above. 

The outline is as follows. The mean- field 4 model 
and its basic properties will be discussed in the Sec. [H] 
In Sec. IIIII we will investigate the topological properties 
of the model. In particular, we will prove that the claim 
that topological and thermodynamic transition points co- 
incide is not correct, in general. The final section contains 
discussion where we explain the origin of this discrepancy. 
Some more technical details are given in the Appendix. 



II. MEAN-FIELD 4 MODEL 

Let x n be a scalar displacement of a particle from a 
lattice site n. We consider the following potential energy 



N 

n—1 ~ 




(8) 



depending on the A^-particle configuration x — 
(xi, . . . , xn). Vq(x,H) is an asymmetric on-site poten- 
tial: 



V {x,H) 



-xH 



1 



1 



(9) 



that becomes symmetric for H = 0. The final term in Eq. 
© represents the harmonic interaction between all parti- 
cles with a coupling parameter [i > 0. The reader should 
note that the potential energy and the displacement can 
always be scaled such that Vq(x, H) has the cc-dependence 
given in Eq. ©. This type of models was used to de- 
scribe structural phase transitions [23j . In contrast to the 
mean- field models studied in Refs. |8l ITlL IL2I Il3] | there is 
a nontrivial coupling constant /i, which can not be put 
to one by an appropriate scaling of the temperature. 

Some thermodynamic properties of the model de- 
scribed by Eqs. ||B} and as well as some features 
of its PES were already investigated |2^|. Let us recall 
these results and start with the thermodynamic behav- 
ior. Due to the infinite-range interaction the mean-field 
approximation becomes exact for N — *■ 00. This leads 
to the self-consistency equation for the order parameter 
(x) = (x n )(T,H): 



dxx exp 



Vq(x, H) - fx (x) x 



T 



(10) 
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where Z is the corresponding partition function and j3 = 
1/T. Of course, a phase transition (second order) occurs 
at some T c for H — 0, only. T c follows from: 



jZ o dxx 2 exp[-V {x,0)/T c ] 
1 n o dxexp[-V (a;,0)/T c ] ' 

For < [i -C 1 one finds: 

T c (/x)= M + (9( M 2 ) 



(11) 



(12) 



which yields for the average potential energy per particle 
v(T)= lim N- 1 (V (x)) (T, iJ = 0) at T c 



v c = v{T c ) = --(l-2fi)+0(ii 2 



(13) 



which is always larger than the minimum value —1/4 of 
V (x,0). 

Let us now turn to the stationary points as discussed 
in Ref. El. With the "internal" field 



N 



N 



the solution of dV/dx n reduces to that of 
x 3 - x - H cS = 

with the effective field [2^] 

H c ft = H + flint- 



(14) 

(15) 
(16) 



For |ff eff | < H c = 2/(3^3) there are three real roots of 
Eq. l(TK|l which will be denoted by x a (H c s), a = +, 0, — . 
It is x + > xq > X- . A stationary point of V is charac- 
terized by N/y , the number of x n in x which are equal to 
Za(H c ff). Permuting the particle indices yields station- 
ary points with the same potential energy. Since {x n } 
are displacements and not positions of particles in a liq- 
uid, these permutations should be counted as different 
stationary points. Hence there are: 



P(N+,N ) 



N\ 



N+\N \NJ. 



(17) 



stationary points of class (iV + , iVo, iV_), where 

^N a = N, (18) 

cr 

i.e., iV_ — N — N + — No. The characterization of all 
stationary points by iV_|_, No and iV_ or equivalently by 



= ^ V. =1 

a AT ' 2-^i 1 



N 



(19) 



proves (see, e.g., Sec. 3) to be extremely useful. Hav- 
ing specified a a we can determine i/eff (ct+, a o) from 
Eqs. (O and {TBJ}: 



H cff = H + fi^2 at a x a (H eS ). 



(20) 
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FIG. 1: a - (^-dependence of v(a+,ao) for n = 0.2 and 
H = 0; b - Same for H = 0.1. Dotted lines locate the maxima. 
Note the concavity of v(a+,cto) in q+ for H > and its 
symmetry with respect to the maximum position a™ ax = (1 — 
a )/2 for H = 0. 



[We omit the arguments H and \i for brevity and we 
take into account that a_ can be expressed by a+,ao, 
due to Eq. I|19f> . Also we will mostly drop the arguments 
(cv + ,ao) of H c ff]. Finally the stationary points of class 
(N + ,N , N_) are given by: 



x(H c ff) = ( x + (H c fi), ■ ■ ■ ; xo(ffcff), ■ ■ ■ ; x-{H c s), ■ • ■ ). 

N + times No times iV_ times 

and its permutations. Similarly for the potential energy 
per particle u(a+,ao) follows: 

v(a+, a ) = ^2 a a Vo{x a (H cff ),H) - — [H e g - H} 2 . 

(21) 

Fig.^presents u(aq_, ao) as a function of a+ for different 
values of the parameters. 

Eq. I|21|) holds for \H c g\ < H c . One can easily prove 
that the latter is guaranteed for arbitrary {a a } with 



4 



1 if 



-H c (l - 3/x) < if < (1 - 3p)H c . 



(22) 



This implies the necessary condition /i < 1/3 (cf. also 
Ref. |24|1. fx = 1/3 is a critical value for the coupling 
parameter at which the PES changes qualitatively. If 
inequality Eq. l|22|l is violated, the stationary points are 
not uniquely classified by {a a } 26]. Therefore we will 
restrict ourselves to /i < 1/3, in the following. What 
remains to be done is to determine the saddle index n s 
of these stationary points. The answer is simple [24| . 
because 



a . 



(23) 



This has been proven for H = in Ref. [24| by determi- 
nation of the number N s of negative eigenvalues of the 
Hessian at stationary configurations with Ao fixed. The 
result is that N s = No, which implies Eq. I|23|l . Validity 
of Eq. 1)23(1 can also be seen as follows. If /i = and 
H > 0, the stationary points of V are given by those of 
Vb and the sign of the eigenvalues of the Hessian equals 
the sign of V . Since Ao is the number of particles with 
position xq(H), which is on the concave part of Vq, i.e. 
V (xo(H)) < 0, it is No = N s . N s is a topological in- 
variant for < fi < 1/3. Consequently Eq. (|23|l remains 
true for all /i < 1/3 and all H obeying Eq. J23 ■ 



III. STATISTICS OF STATIONARY POINTS 

A. Saddle index distribution 

We have shown in Sec. [H] that all stationary states 
of the mean-field <fi 4 model are characterized by {a CT }, 
provided inequality Eq. I(22|l holds, ao = n s is the sad- 
dle index. Consequently there are in total 3 N stationary 
points, from which 2^ are local minima. The saddle in- 
dex distribution follows from Eq. i|17|) : 



N-N s 



p(n s ) = 3- N J2 p ( N +' N s 



-NnN — N s 



N\ 



N+=0 



N S \(N-N S )\ 



(24) 



Using the Stirling formula for N 1 one obtains 

p(n s ) ~ exp [Ns(n s )] , 
where the configurational entropy is given by 

1 - n. 



s(n s ) = - 



n s \nn s + (1 — n s ) In ■ 



(25) 



The maximum of p(n s ) is at n™ ax = 1/3 which is obvious 
since P(N + ,No) has a maximum at N + = No = iV_ = 
N/3. It is interesting that this result based on the double- 
well character of the local potential, coincides with the 
numerical finding for binary Lennard- Jones clusters |17| . 



We stress that the validity of n™ ax = 1/3 is more general. 
Suppose there are no interactions. Then the iV-particle 
problem separates into that of N independent particles 
in double wells, for which P(N + , No) is still given by 
Eq. (|17f) . Turning on an arbitrary interaction will not 
destroy the one-to-one correspondence between station- 
ary points and (JV+, -/V_, iVo), up to a critical interaction 
strength. At this critical strength, e.g., exponentially 
many metastable configurations may become unstable. 
Accordingly, n™ ax = 1/3 holds up to that critical cou- 
pling, i.e., it is a topological invariant. 



B. Calculation of the density of states p(e,n s ) 

The density of stationary points of a PES with energy 
e and saddle index n s that was mentioned in the Intro- 
duction is defined by 

JV-JVo 

p(e,n s )= P(N+,No)5(v(a + ,ao)-e), (26) 

N + =0 

where P(N + ,N ) is given by Eq. l(T7|l and 

N = a N = n s N, N + = a + N. 

The density p(e) of stationary points with energy e fol- 
lows from 



N 



'(e) = ^2 P( e > n *)> 



(27) 



Af o =0 



Neglecting the irrelevant prefactor one can immediately 
write (ao = n s ) 



p{e,n s ) - P(N + (e,a ),N ), 



(28) 



where N + (e, ao) — a+(e, ao)N follows from the solution 
of the equation 



v(a + ,a ) = e. 



(29) 



This equation has two solutions a^(ao, e) which are de- 
rived and discussed in Appendix B. 

With the use of the Stirling formula P(N + ,N ) sim- 
plifies to P(N+, No) ~ exp [Ns (e, ao)] , thus one obtains 
Eq. @ with the configurational entropy given by 

s (e, ao) = —ao In a — a+ In a^ — a^ In a^ , (30) 

where a^ 1 = 1 — ao — a+ and we have taken the (+) 



branch of Eq. (|B1|) that makes the dominant contribution 
into P{N + ,N ) for H > 0. For H = one can also 
restrict oneself to a^ , due to the symmetry. Note that 
p(e, n s ) is nonzero and given by Eq. @ only in the energy 
window 



Wmin(ao) < e < Umax(ao), 



(31) 



5 




where v min (a ) = v{l-a ,a ) (see Fig. [TJ and u max (a ) 
is given by Eq. (|A7|I . otherwise p(e,n s ) = 0. Alterna- 
tively one can say that Eq. is valid in the window of 
saddle indices 



(min) 



/ \ (max) / \ 

(e) < a < a K '(e), 



where the boundary values satisfy v m i n (ctQ 
w ma x((to ml °') = e - From Eq. I|A7(1 one finds 

4 min) (e) = l + 4e-2H 2 /», 



(max) 



(32) 



e and 



(33) 



whereas (e) can be found with the help of Eq. I|21|) 

or, approximately, with the help of Eq. (|A5|I . The de- 
pendence of s(e, n s ) on n s is shown in Fig. Elfor zero and 
nonzero field. 

Let us discuss the main features of s(e, n s ) presented 
in Fig. |21 For a more detailed analytical discussion of 
s(e,n s ) the interested reader is referred to Appendix B. 
We begin with H = (see Fig. |2K)- The maximum 
of s(e,n s ) with respect to n s is denoted by n s (e). Be- 
cause of the relation between p(e,n s ) and s(e,n s ) given 
by Eq. © it is obvious that for N — > oo the maximum 
position n s (e) is identical to the averaged saddle index 



n s (e) given by Eq. Q). In Appendix B the existence of a 
critical energy e c (/z) is proven. s(e,n s ) as function of n s 
has a maximum within the domain of n s for e < e c ([i) and 
a maximum at the left border of its domain for e > e c (/x). 
This implies that the slope ds(e,n s )/dn s at n s (e) is con- 
tinuous in e, but not differentiable, i.e., the "curvature" 
d 2 s(e 1 n s )/dn 2 is discontinuous in e at e = e c (p). This 
is the origin of the topological singularity, discussed be- 
low. Fig. 3a presents n s (e) and reveals the singularity 
at e c (/i). Note that n s {e) contains a branch that is inde- 
pendent of the interaction [see Eqs. (|B6|I and l|B7|) ]. For 
e very close to the ground state energy vq(/j,) [cf. i|A6(l ] 
one obtains the power-law behavior (see Appendix B): 



n s (e) 



Ml*)} 



<5(M) 



(34) 



with 5(p) > 1, if n is small enough. Note that n s (e) — > 
for e — > v (p) ■ The averaged saddle index n s (e) takes the 
critical value ni c \fi) — n s (e c (p) shown in Fig. 0] This 



figure also includes the asymptotic result of Eq. l)B5jl for 
fi 0. 

Instead of fixing e, one can also determine the maxi- 
mum position e(n s ) of s(e,n s ) for given n s . e(n s ) is the 



6 




n s (e) 




FIG. 4: ^-dependence of the critical average saddle index 
rif\n) for < /i < 1/3 (solid line). The dashed line is the 
asymptotic result Eq. <B5t for fj, <C 1. 



averaged saddle energy [cf. Eq. J5J] as function of n s . It is 
easy to get e(n s ), since the number of stationary configu- 



rations is maximal for ojj 



(l-n s )/2. This yields 



for the effective field H c s(a+ = (1 — n s )/2, ao = n s ) = 
which implies x±{a+ = (1 — n s )/2,ao = n s ) = ±1 and 
xo(a+ = (1 — ri s )/2, ao = n s ) = and this in turn leads 
to 



e(n s ) = v(a + = (1 - n s )/2, a = n s ) = -(1 - n s )/4. 

The reader should note that (i) the inverse function 
n s {e) (see Fig. of e(n s ) turns to zero at e = —1/4 
which equals the lowest energy of the on-site potential, 
but is above the ground state energy vo(fi) and (ii) e(n s ) 
is not the inverse of n s (e). This difference is related to 
the fact that the maximum of p(e, n s ) with respect to n s 
for fixed e is not generally related to its maximum with 
respect to e for fixed n s . 

Now let us take H ^ (see Fig. |3Jd). As discussed in 
Appendix B, s(e, n s ) has always a maximum at n s (e) as 
function of n s within its domain. n s (e) is shown in Fig. 
|3Jd. The nonsingular e-dependence for H ^ can clearly 
be seen. The average energy e{n s ) or its inverse n s (e) 
can be derived analytically (see Appendix B). The result 
for n s (e) is shown in Fig. |SJ 

Having determined n s (e) for H = and H ^ 0, one 
can now calculate the energy-dependent configurational 
entropy that follows from Eq. I|27|l . in the limit N — > oo 



s( e ) = s(e,n s (e)). 



(36) 



The result is shown in Fig. |SJ As can be seen from Fig. 
|BJ), s(e) has a discontinuous second derivative for H = 
at e = e c (/i). In Fig. these points are marked by 
circles. The high-energy branch of s(e) has the form of 
Eq. H25[) with n s (e) =>• 1 + 4e. It attains a maximum for 
n s (e) = 1/3 that implies e = —1/6, independently of the 
interaction. 



FIG. 5: Saddle index n 3 (e) versus averaged energy e for 
H = 0.2 and H = and 0.15. 



C. Euler characteristic x( e ) 

It turns out that for the model under consideration in 
the limit N — > oo the Euler characteristic of Eq. 
satisfies 

lx(e)|~p(e), iln|*(e)|=*(e). (37) 

Calculation of x(e) for large N is similar to that of s(e), 
as suggested by the similar form of Eqs. J3J) and The 
only difference is that M(e, 7V S ) used in the calculation of 
X(e) is, unlike p(e, n s ), nonzero for v max (ao) < e. In this 
case it is independent of the energy and has the form 



M(N S ) = 



■yN-N. 



N S \(N-N S )\ 



(38) 



similar to Eq. (|24|l . Hence there are two different con- 
tributions into x(e). It can be easily shown that the con- 
tribution from the range v m in(Q!o) < e < v max (ao) for 
N — > oo coincides with that of p(e) ~ exp [iVs(e)] that 
was studied above up to an irrelevant prefactor, in spite 
of the sign alternation in Eq. QJ. In contrast, the contri- 
bution from the range of N s determined by w max (ao) < e 
is dominated noi by the maximum of M(N S ) on iV, but 
by N s on the boundary of its interval, i.e., by N s sat- 
isfying v max (N s /N) — e. One can easily see that this 
contribution to x( e ) never exceeds that from the range 
^min (ao) < e < -L> max (a ), thus one obtains |x(e)| — p(e) 
for — > oo. 

The reason for such a peculiar behavior of the contribu- 
tion from the region f m ax(ao) < e is the sign alternation 
in x(e) plus the specific form of M(N S ), For instance, as 
all stationary points are below the level e = 0, one finds 
x(e) for e > 0, by just summing over all N s : 



X(e > 0) 



N 

E 

N a =0 



(-1)">M(N S ) = 1. 
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FIG. 6: a - Energy dependence of the configurational entropy 
s(e) for different interaction strength fj,, with and without 
field. Dashed line is the solution for /i = and H = and 
the circle indicates the location of e c (/x); 6 - Derivative of the 
configurational entropy showing a transition at e = e c (fj,) for 
H = 0. 



This result is exact and it has a transparent topological 
meaning. Replacing the sum by the maximal summand 
value M(N/3) would be an error. Even simplifying Eq. 
(I38|l with the help of the Stirling formula for N ^> 1 in the 
sum would lead to an exponentially large result instead 
of 1. Therefore, one should be cautious in applying the 
saddle point method to the r.h.s. of Eq. 0J. 



IV. DISCUSSION 

We have investigated the statistics of stationary points 
and topological properties of the analytically tractable 
potential energy surface of a (ft 4 model in a symmetry- 
breaking field H with interaction of all pairs of particles 
with the same strength /i. For this model the mean- 
field approximation becomes exact in the thermodynamic 
limit N — > oo. For H = there is a second-order phase 
transition at the critical temperature T c ([m) that is ana- 
lytical in [i. 



We have shown that the distribution of the saddle in- 
dices p(n s ), where n s — N s /N and N s is the number of 
unstable directions at a stationary point has a maximum 
at n s — rt™ ax = 1/3. Interestingly this value is consistent 
with that found for small binary Lennard-Jones clusters 
[l7j . Whether or not this is an accident is not clear. Our 
result originates from the fact that all stationary points 
can be labelled by symbolic sequences (<7i, • • • , cr n ) with 
<r„ = +,0,—. The low-temperature anomalies of struc- 
tural glasses are usually explained by the existence of 
two-level systems arising from an ensemble of double- well 
potentials. As the smallest "unit" of a PES of a classi- 
cal ./V-particle system, one may choose the local minima 
including their basins of attraction. But such a choice 
does not fully encompass the saddle. Taking the next 
larger unit, a pair of local minima and their common 
saddle, one arrives at a double-well characterization of 
the PES. This could explain why n™ ax = 1/3 for small 
clusters and for liquids. We have also argued that the 
value n™ ax = 1/3 is a topological invariant for an entire 
family of (ft 4 models. In any case, it would be important 
to determine n™ ax for other e.g., liquid- like models and 
to check whether it equals again n™ ax = 1/3. 

For our model the absolute value of the Euler char- 
acteristic x( e ) is essentially the same as the density of 
stationary points p(e) in the limit N — > oo, see Eq. I|37ll . 
It would be interesting to investigate the generality of 
this result. 

For H = we have found a singularity in p(e) and thus 
in |x(e)| at the energy e c (/z) given by Eq. i|B5|) . e c (/i) 
is nonanalytic in fi. At e = e c (ji) the second derivative 
d 2 In \\\ (e)/de 2 is discontinuous, as found for the models 
studied in Refs. |lfJ . lllLll2lll3| . In agreement with these 
papers, we also have found that the topological singu- 
larity disappears for nonzero field, as well as the ther- 
modynamic singularity. In this respect, we would like to 
mention a recent publication Ref. [27j , where it is proven 
that a topological singularity is a necessary condition for 
a thermodynamical transition to take plac e. However, at 
variance with earlier work 0, ITTl IT^L Il3| , this singular- 
ity is not related to the thermodynamic singularity of our 
model as e c (/i) does not coincide with v c (/j,), the average 
potential energy at temperarure T c (fi). 

The reasons for this discrepancy can be made clear 
with the help of the following argument. Let us 
smoothly modify the local potential Vq(x) in the inter- 
vals [-oo,-l- e(fj.)], [-l + s(n),-e(fx)], [e(/i), 1 - e(ji)] 
and [1 + e(/x),oo] for given < e(/i) < 1. If p, is small 
enough, which implies that the internal field is small 
enough, then the three roots of Eq. H15[) are within the 
intervals [a — e(/j,),a + e(/i)], er = +,0,— , in which the 
potential has not been modified. Accordingly, the roots 
and therefore the stationary points and their energies are 
the same. This implies that e c (fi) is the same. However, 
since the calculation of T c (p) [cf. Eq. JTTJ] involves Vq(x) 
for all x, the critical temperature will be different for the 
modified on-site potential. 

The discrepancy between the topological and thermo- 
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dynamic singularities can also be traced back to an un- 
justified comparison of a continuous model (thermody- 
namics) and a discrete model (topology) . More logically, 
the energy of all the stationary points can be represented 
by an Ising-like Hamiltonian H({(Ji}) (for /i < 1/3) with 
<Ti = +, 0, — . The corresponding canonical partition func- 
tion Z(T) = Tr exp [— 7i({(Tj})] can be calculated from 
the density of states p(e) as 



Z(T) = J dep(e)exp [-Ne/T]. 



(39) 



Evidently a singularity of Z(T) at the corresponding 
transition temperature T' c results from the underlying 
singularity of p(e) at e c . Obviously in this case v' c (n) = 
(H)(T' C )/N = e c (/i) is fulfilled. But the thermodynamics 
of this discrete model does not coincide with that of the 
original continuous model, in particular, T' c ^ T c . 

The idea of an at least qualitative relationship be- 
tween the thermodynamic singularity and the topolog- 
ical singularity is supported by the following observa- 
tion. At the thermodynamic transition point there ap- 
pears a spontaneous breaking of the left-right symme- 
try for the displacements, which is equivalent to emerg- 
ing of a nonzero temperature-dependent internal field for 
T < T c (/i). On the other hand, the stationary configu- 
rations with e > e c (fi) and with maximum weight corre- 
spond to a + = a_. This implies that the effective field 
defined by Eqs. lfH|) and l(TB|l satisfies if ff = 0. How- 
ever for e < e c (//) it is a+ ^ a_ and hence H c g ^ 0. 
Therefore a spontaneous symmetry breaking occurs at 
both singularities. 

The energy- or temperature-dependence of the aver- 
aged saddle index n s seems to play a role for the dy- 
namical features of supercooled liquids. For the present 
model we have found that n s (e) vanishes at the ground 
state energy vq, only. This is consistent with recent re- 
sults on n s (T) showing that n s — at T = 0, only [l3j . 
Taking the analogy to mean-field like spin glass mod- 
els, this would imply that no dynamical transition (or 
crossover) could occur at finite temperatures 22] . The 
averaged energy e(n s ) as function of n s is not the in- 
verse function of n s (e). In particular n s (e) vanishes at 
e*(H) > vq. Whether dynamics changes qualitatively at 
e*(H) or any other characteristic temperature would be 
interesting to study. 

Finally we would like to mention that after submission 
of this paper we learned about a similar study of the same 
model 28], where the authors find, among others, the 
same type of discrepancy between the thermodynamical 
and topological singularities. 



APPENDIX A: PROPERTIES OF v(a+,a ) 

In this Appendix we investigate the properties of 
v(a + ,a ) defined by Eq. J2U- For the discussions in 
Sec. 3, the first and second derivative of v(pt.+ , ao) with 



respect to a+ will be useful. A straightforward calcula- 
tion making use of Eqs. I|15|). l(Pp|l and (|2U|) yields: 



dv(a + ,a ) 
9a + 



V Q (x + (H eff ), H cS ) - Vo(x + (H cS ), i^) 



(Al) 

Since H c s = implies x± = ±1 [see Eq. Q15[l] it follows 
that dv(a+,ao)/da + = for H e g = 0. That is, the 
maximum of v(a + ,ao) with respect to a+ corresponds 
to ff c ff = 0. For the second derivative we get: 

d 2 v(a + , a Q ) 



da 



2 = -[x+{H eS ) - X-(H eS )] 

dH cS (a + ,a ) 
da + 



(A2) 



with 



dH cS (a + ,a ) _ x + (H cS ) - x-(H cS ) 

— aZ ■ \ A6 ) 



da+ 



1 ME 3x-i(Hln)-l 



Since x+ — x_ > 0, the "curvature" d 2 v/da 2 i _ is negative, 
i.e., v(a+, «o) is concave in a+, provided dH c s/du + > 0. 
This is true if fx is small enough, as can be seen from 
Eq. fftSl . 

The function v(a+,ao) can be computed analytically 
near its maximum in a + by using the expansions 

x-(H eS ) = -l + H eS /2 
x (H cS ) = -H eB 

x+(H cf{ ) = l + H cS /2 (A4) 
near H c h = 0. The result is a parabola in a + : 



1 H 2 

v(a + ,a ) = --(1 - a ) + — 

1 [H + m(2q+ + a - l)f 
2/x l-/x(l-3ao)/2 



(A5) 



Analysis shows that corrections to this formula in the 
whole region < a + < 1 — ao are of order /i 3 , /j, 2 H, 
and /J,H 2 , i.e., Eq. I|A5[) is a very good approximation for 
not too large [i. For instance, the ground-state energy 
following from Eq. i|A5|l in the case H = 



W = 0(0,0) = -i(l +T -?^ 

is in accord with the exact two-fold degenerate ground- 
state energy, 

u (aO = "(1, 0) = "(0- 0) - ~(1 + y) 2 . (A6) 

up to the terms /i 2 , and the relative error is only 0.0125 
for fi = 1/3. The exact value of the field-dependent max- 
imal potential energy that also follows from Eq. (|A5|) has 
the form 



1 H 2 
fmax(ao) = (1 ~ ao) + y ■ 



(A7) 
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The corresponding exact value of a + is 

(max), n 1 - OLQ H 

Eqs. l|A7jl and (|A8j) are thus valid for \H\ < fi(l - a ). 



Its solution a.Q (fi) = nf' (/i) that is plotted in Fig. 0]and 
the corresponding energy e^ c \fi) are critical parameters 
that define the boundary between different regimes. In 
particular, a c) (l/3) ~ 0.1774 and e( c )(l/3) ~ -0.2056. 
For fi -c 1 one can solve Eq. 1B4|) analytically: 



APPENDIX B: CONFIGURATION AL ENTROPY 

In this Appendix we present analytical results for the 
configurational entropy s(e,n s ). In particular, for H = 
the existence of a critical energy e c (/x) will be proven, at 
which singular energy dependence occurs. 

First we have to investigate v(a+, ao). It is convenient 
to use the approximate form of u(a+,ag) given by Eq. 
(|A5I) . Then Eq. i)29|) becomes quadratic, and one obtains 
two solutions 



„(±) 



'(«o) 



>maxM-e] [1 - m(1 - 3oo)/2] 



2/y 



,(B1) 



where a^ 1 (ao) and v 



max (ao) are given by Eqs. (SHI 
Note that this result becomes 



and (S3 ) respectively, 
exact for all /i with < // < 1/3 provided that e is close 

tO U m ax(ao)- 

The density of stationary points p(e, n s ) can be used to 
calculate the statistics of saddle indices, see Eqs. Q-©. 
The value of n s (e) that maximizes the entropy s (e,ao) 
on ao is the solution of the equation 



Oa 



(+) 



da 



■In 



1 — ao — a 



(+) 



In 



1 — ao — a 
ao 



(+) 



= (B2) 



for ao, if this solution exists in the interval of Eq. (|32ll . 
Otherwise the maximum of s (e, ao) is attained at the left 
boundary of the ao-window, ao = ag (e) (see Fig. |2Ji)- 
The latter solution exists only for H = 0. It corresponds 
to the maxima of the curves w(a+, ao) in Fig. {Tfr, i.e., to 



a 



(+) 



(max) 



(«o) 



1 - a 



(B3) 



One can see from Eq. (|B1(I 
e -> Umax (a ) , i-e., for a 



that da 

(mill 



(+) 



/<9ao diverges for 

i.e., tor ao — » ag '(e). One can check 
that in the case H = this divergence is compensated for 
by the log factor that tends to zero. For H ^ there is 
no such compensation, and the lhs of Eq. I|B2(1 diverges 

at ag — > a mm ' ) (e), thus Eq. 1B2I) always has a solution, 
see Fig. 

Let us consider the case H — and find the condition 



that the solution of Eq. I|B2|) 
Simplifying this equation for e 
the transcedental equation 



is just ag = a 



(min) 



(e). 



(ag) one obtains 



1 - /i(l - 3a )/2 1 - a 



4/x(l - a ) 



In- 



2a 



(B4) 



e c (M) = -- 



'(M) 



2 

1 1 

-4 + 8 



4/i 



exp 



1 

4^ 



(B5) 



Note that e^ c \fj,) is always above —1/4, the ground-state 
energy without interaction. Now one can write down the 
combined expression for n s (e) in the case H = 0: 



.(e) 



n™ ax (e), 



wo(a*) < e < e c (/i) 
e c (/x) < e < 
< e. 



(B6) 



Here Vo(/i) is the ground-state energy, nf(e) is equal to 
a that solves Eq. (|B2() . and 



nr x (e) = 1 + 4e. 



(B7) 



Note that the high-energy branch n™ ax (e) is independent 
of fj, and is thus the same as for a system of noninteract- 
ing particles. This contribution is due to the maxima 
of v(a + ,ao) [see Eq. |B3jl] whereas nf(e) is the contri- 
bution from the energy levels below this maximum. One 
can show that n s (e) — * for e — * vq(^i). This dependence 
has the form given by Eq. I|34|) . However this dependence 
is only realized for the energies very close to the ground 
state vo(fj,). The function n s (e) is shown in Fig. [31 for 
H = and if ^ 0. It has a discontinuous derivative at 
e = e c (/i) in zero field. This discontinuity disappears for 

One can also use p(e,n s ) to calculate the average en- 
ergy e(n s ) for a given n s . From the second of Eqs. © 
one obtains 



1 - a 



a 



(+) 



„(+) 



In 



that is simpler than Eq. (jB2j|. From its solution 



e(n s ) 



H 



1 — 3n» 



4 l-/i(l- 3n s )/2 



(B8) 



one obtains n s (e) shown in Fig. [SJthat is almost linear in 
e in the considered range of parameters, Eq. (|22|l . Note 
that, in contrast to n s {e), the quantity n s {e) found from 
Eq. (|B8|I turns to zero not at the ground-state energy 
but at e = -1/4 - H 2 / [2(2 - fi)]. 
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